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The A + B — » C reaction-diffusion process is studied in 
a system where the reagents are separated by a semiperme- 
able wall. We use reaction-diffusion equations to describe the 
process and to derive a scaling description for the long-time 
behavior of the reaction front. Furthermore, we show that 
a critical localization-delocalization transition takes place as 
a control parameter which depends on the initial densities 
and on the diffusion constants is varied. The transition is be- 
tween a reaction front of finite width that is localized at the 
wall and a front which is detached and moves away from the 
wall. At the critical point, the reaction front remains at the 
wall but its width diverges with time (as i 1//6 in mean-field 
approximation) . 

PACS numbers: 82.20Wt, 82.20Db, 82.20Mj, 66.30.Ny 

I. INTRODUCTION. 

Reaction fronts formed in diffusion- limited A + B — > C 
type reactions have been investigated intensively in re- 
cents years f^-Q. The motivation comes partly from 
the realization that moving reaction fronts play an im- 
portant role in a great variety of physical and chemical 
phenomena which display pattern formation |2^] . An- 
other reason for the interest is the simplicity of the prob- 
lem which allows the application of different theoretical 
approaches. Indeed, the front properties have been stud- 
ied in detail by using mean-field and scaling theories || , 
dynamical renormalization group p[ , numerical simula- 
tions [Q) and in some cases exact analytical predictions 
have also been made Jl3[ . 

In most of the cases studied previously, the reaction 
front is formed after the spatially separated components 
A and B come into contact. For example, in a typical 
experiment aimed at producing Liesegang bands ]25j ] , one 
has a vertical tube of gel soaked with component B, and, 
at time t = 0, a liquid containing the reagent A is poured 
over the gel (in order to eliminate convection effects, the 
liquid is sometimes replaced with another gel containing 
A). The theoretical equivalent of this situation is that 
the reagents are separated by a wall which is removed at 
t = and then the reaction-diffusion process begins. 



One can imagine, however, that there are situations 
when the wall between the reagents is present at all times, 
and this wall is semipermeable allowing only one of the 
reagents to pass through. It may happen, for example, 
in the above discussed setup that B is not soluble in 
the liquid containing A which is effectively equivalent to 
the presence of a semipermeable wall. More importantly, 
chemical reactions in biological systems take usually 
place in strongly inhomogeneous media with semiperme- 
able walls present [p6|-p8[. Thus, we believe it is impor- 
tant (hence the aim of this paper) to consider the forma- 
tion of reaction fronts in systems with initial separated 
species when the wall separating the two species is not 
eliminated at t = but is replaced by a semipermeable 
wall which allows only one of the reagents (A) to diffuse 
across. 

Using a mean-field description of the above process, we 
find that the control parameter in this system is given by 

r = l 7=- (1) 

a VD a 

where eto and bo are the initial particles densities while D a 
and Db are the diffusion coefficients. We show that, de- 
pending on the sign of r, three distinct types of behavior 
occur. When r > 0, the A particles invade the B phase. 
The reaction front moves away from the semipermeable 
wall with the distance from the wall increasing as \ft and 
the wall is irrelevant in the long time regime. Thus, one 
recovers the predictions (e.g. the width of the reaction 
zone scales as id ~ i 1 / 6 ) made with no semipermeable 
wall present p]. In the opposite case, r < 0, the wall 
prevents the B particles from invading the A region and, 
accordingly, the reaction front becomes localized (with fi- 
nite width) at the semipermeable wall. It turns out that 
the dividing point between the r > and r < cases 
is a critical point in the sense that the width of the re- 
action zone diverges at r — 0. We have thus found a 
critical localization-delocalization transition from a reac- 
tion front localized at the wall to a front detached and 
moving away from the wall. 

The above results will be derived and discussed first by 
defining the model (dynamical equations and the bound- 
ary conditions) in Sec. II. Then, the different regimes are 
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analyzed (Sec. Ill) both analytically and numerically at 
the mean-field level with comments on the role of the 
fluctuations. Concluding remarks are given in Sec. IV. 



II. THE MODEL 

The basic notions about reaction zones have been in- 
troduced for the A + B — ► C process |j] and, in or- 
der to keep the discussion transparent, we shall also 
consider this case. More complicated reaction schemes 
uaA + vbB — > C can be treated along the same line with 
the same general picture arising. 

We shall assume that the transport kinetics of the 
reagents is dominated by diffusion and that the reaction 
kinetics is of second order. Thus, at a mean-field level, 
the mathematical description of the process is given in 
terms of reaction-diffusion equations 



and the width of the reaction zone 



dta = D a \7 2 a — kab 
d t b = D b V 2 b - kab 



(2) 
(3) 



where a and b are the densities of the reagents A and B, 
respectively, D a and Db are the corresponding diffusion 
constants, and the reaction-rate parameter is k. Note 
that there is a conservation law in this system. Since 
the A and B particles react in pairs the difference in 
their numbers is conserved. In terms of the densities this 
means that the spatial integral of a — b is constant unless 
there are particle sources at the boundaries. 

The semipermeable membrane is located at the (x — 
0, y, z) plane. Initially, all B particles are on the right 
hand side of this membrane (x > 0) and, since the mem- 
brane is impenetrable for them, they remain on that side 
for all times. In terms of the particle density b this means 
that the solution of (j|) and ||) must satisfy the following 
conditions 



b(x < 0,t) = 



db(x,t) 



dx 



= 



(4) 



2 = 0+ 



The motion of the A particles is not influenced by the 
membrane and, initially, they are on the left side of it. 
Furthermore, the initial densities are assumed to be con- 
stant i.e. a(x, 0) = ao and b(x,0) = for x < while 
a(x, 0) = and b(x, 0) = b for x > 0. With this choice 
of initial state, the solution of (j|) and (||) depends only 
on the x spatial coordinate and the system effectively 
becomes one-dimensional. 

Our aim will be to calculate the production rate of C 
particles defined by 



R(x,t) = ka(x,t)b(x,t) 



(5) 



and investigate the time-evolution of its spatial structure 
with emphasis on the center 



/oo />oo 
xR(x,t)dx/ / R(x,t)dx 
-oo J —oo 



(0) 



w(t) 



(x — xf) 2 R(x,t)dx/ R(x,t)dx 



1/2 



(7) 



Both Xf and w are the easily measurable quantities in 
experiments and simulations. 



III. SCALING PROPERTIES OF THE FRONT 




For a system without the membrane, it is known [p||l7|] 
that the reaction front will move to the right (A invades 
B) or to the left (B invades A) depending on the relative 
magnitude of quasistationary diffusive currents (J A ~ 
D a ao/\/D a t and J B ~ Dj,bo / \fD^t) , i.e. depending on 
the sign of the control parameter r: 



(8) 



For r = r c = 0, the front is stationary in the sense that 
although t) remains time-dependent for large times, 
the center of the reaction zone does not move and Xf(t—> 
oo) approaches a finite constant. 

One expects that the direction of invasion plays an 
important role in the presence of the membrane as well 
and, accordingly, we shall analyze the r > 0, r < 0, and 
r = cases separately. 



A. r > 0: Invasion of the free (A) reagents — 
delocalized front 

For r > 0, the diffusive current of A particles over- 
whelms the corresponding current of B particles and thus 
the reaction front moves to the right. After a while, the 
B particles disappear from the neighborhood of the mem- 
brane and thus the membrane does not play a role any- 
more. Consequently, the reaction front leaves the mem- 
brane (Fig.|l]) and all the results about the long-time scal- 
ing form of the reaction front obtained previously ap- 
ply H namely 



R(x,t) - r f3 F 



x f 



t a 



(9) 



where the position of the center of the front, Xf, scales 
with time as Xf ~ y/i, the width of the reaction front is 
proportional to t a with a = 1/6 , the scaling exponent 
of the production rate of C at x = Xf is (3 — 2/3, and 
the scaling function, F(z), is a fast decreasing function 
for z — > ±oo. 
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FIG. 1. Production rate of C near the semipermeable mem- 
brane (localized at x = 0) for three different values of the 
control parameter r. The width of the localized front is sta- 
tionary while the width of the critical front increases with time 
as i 1 / 6 . The distance between the delocalized reaction front 
and the membrane increases as t 1 ^ 2 while its width diverges 
as t 1 ' 6 . 

We can call this front delocalized since both the center 
and the width of the front diverge in the long-time limit. 

In closing this subsection, we note that the above 
results are modified by fluctuations in low dimensions 
(d < 2), as discussed in several works on A + B — > C 
reactions without the presence of a membrane |],[ll| . 

B. r < 0: Invasion of the blocked (B) reagents — 
localized front 



the reaction zone. The second condition in ( [l0| ) is more 
complicated. It follows from the assumption that the 
penetration length, £, is finite combined with the fact 
that the diffusion current approaches zero at large times, 
[Jdiff ~ D a (da/dx)\ x= a — ► \j\ft i.e. the derivative 
(da/dx)\ x —o diminishes for t — > oo]. The finiteness of 
£, in turn, follows from the finiteness of 6(0, t) — b* and 
so, finding 6* finite at the end of our calculation provides 
a selfconsistency check of the underlying picture. 
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FIG. 2. Density profile of the reagents for r < as seen 
on a diffusive scale (x ~ t 1 ^ 2 ). Time is measured in units of 
r = 0.1/(fcao) where k is the reaction rate and ao is the initial 
density of A. The unit of length is chosen to be t = \fD a T 
where D a is the diffusion coefficient of A. For the given values 
of diffusion coefficients (D a ,Db) and initial densities (ao,&o), 
the large time limit of b at x = is given by b* — 1 — 1/ y/2. 



For r < 0, the B particles would be the invading parti- 
cles but they cannot penetrate past the membrane. Thus, 
one expects that there will be a finite density of B parti- 
cles at x = and, consequently, the A particles can pen- 
etrate into the x > region only up to a finite distance, 
£. In order to make this picture (Fig.||) quantitative, we 
shall first solve the problem on the diffusive lengthscale 
x ~ y/t and then use this solution as the large-argument 
asymptotics of the solution around x = 0. 

Viewing the process on the diffusive lengthscale, the 
reaction zone is reduced to a point (x = 0) and the diffu- 
sion of A and B takes place separately in the x < and 
x > regions. The appropriate boundary conditions are 
as follows: 



cl^x — ► — oo, = clq , a(0, t) = 



b(x — > oo, t) = b 



-D — 

ox 



x=0~ 



dx 



(10) 



x=0+ 

(11) 



The first boundary conditions in ( |To| ) and ( pi]) a re obvi- 
ous while the second boundary condition in (p_l[) is just 
the expression of the equality of the currents entering 



The solution of the diffusion equations with the above 
boundary conditions is given by 



a(x, 

b(x, t)=b* + (b - 6*)Erf (x/ y/4BM) 



(12) 
(13) 





where Erf(x) is the error function [p9| a nd b* = 6(0, t) is 
found from the second condition in (II 1|): 



(14) 



As we can see, 6* is indeed finite for finite r < (6* > 
because it has the meaning of particle density). 

The above results are valid on lengthscale x ~ t 1 / 2 . 
In order to investigate the details of the reaction zone, 
we must consider the x ~ t° region (Fig.^) where we 
should find a solution with large-distance asymptotics 
which matches to the solution found above. 
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FIG. 3. Magnified view of the reaction zone shown in 
Fig. 2. Here the x coordinate is not scaled by t 1 ' 2 . The pene- 
tration length of particles A into the B region is shown by £. 
Note that there is a break in the vertical scale. 

Since we are mainly interested to find the extent of 
the region where the reaction product appears, we should 
find the region of penetration of A particles into the x > 
halfspace. For x <C one can approximate b(x, t) b* 
and then the equation for a becomes linear 



d t a = D a V 2 a - kb*a 



(15) 



This equation is supplemented with the following bound- 
ary conditions 



a{x^oo,t)=0 , — 
ox 



«() 



(16) 



The second condition comes from the fact that the diffu- 
sion current entering the reaction zone at x = must be 
equal to that calculated from the macroscopic (x ~ \/t) 
considerations. 

Due to the slowness of diffusion, a(x, t) changes slowly 
at large times and one can consider quasistatic approxi- 
mation. Looking for a solution of the form 



a(x, t) 



1 



(17) 



one can see that the left hand side of (|T^) is of the order 
i -3 / 2 while the right hand side is proportional to t -1 / 2 
and so, the time derivative can be neglected. The result- 
ing equation for $ can be easily solved and the boundary 
conditions can be satisfied yielding a solution in a scaling 
form: 



a(x,t) 



^(x/^D a t/e 



e « 



(18) 



a ' ' ' y/wDat/? 

where the penetration (or correlation) length is given by 




-1/2 



(19) 



Since b(x,t) ~ b* in the reaction zone, we can obtain 
R(x,t) from (|f): 

R(x,t) = kab^ kab* ~ X>Q ( 20 ) 



= 



x < 



(21) 



Thus the reaction rate goes down with time as 1/y/t 
while the center and the width of the reaction zone re- 
main finite in this scaling limit 



Xf ~ w ~ £ . 



(22) 



One can see from Fig.|| that the scaling form jl^ ) agrees 
with the scaling obtained from the numerical solution of 
the full set of equations @) and (||) . 
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FIG. 4. Scaling of the density of A's in the reaction zone 
shown in Fig. 2. The parameters are the same as in Fig. 2 
except for t and bo which are varied in order to keep t/£ 2 con- 
stant [£ is given by equation (19)]. The numerical solution of 
the full set of reaction-diffusion equations (equations (2) and 
(3)) is compared with the quasi-stationary scaling solution, 
$ a (solid line). Since one has b w 6* in the reaction zone, 
the scaling function of the reaction rate R = kab =m kab* is 
proportional to that of a. 

The phase considered above may be called the phase of 
localized reaction zone. One can observe from (|l9|) , how- 
ever, that £ diverges as we approach the r — point and 
thus the reaction zone becomes delocalized at r = r c = 0. 
Thus r plays the role of the distance from a critical point 
and the exponent we found, £ ~ r~ v ~ r -1 / 2 is obvi- 
ously the mean-field exponent v = 1/2 in accord with 
the neglect of fluctuations in the above description. 



C. r = 0: Localization-delocalization transition — 
critical front 

It follows from the previous subsection that the r = 
case can be considered as a critical point which separates 
the localized and delocalized phases. Thus we expect 
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that a scaling description is valid again at r = r c = 
but, in expressions like (|l8| ), the correlation length must 
be replaced by a time-dependent correlation length which 
scales as a power of time, ~ t a . In order to see that 
this picture is valid, we follow the steps of the previous 
subsection: the problem is first solved on the diffusion 
scale [the solution is actually given by equations (|l2|) and 
( |l3| ) with b* = 0] and then matching solution in the x w 
region is found (Fig.|and|). 
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FIG. 5. Density profile of the reagents at the critical point 
(r = 0) as seen on a diffusive scale (x ~ t 1 / 2 ). Notation is 
explained in caption to Fig. 2. 




FIG. 6. Magnified view of the reaction zone shown in Fig. 5. 
Note that here x is not scaled by i 1 / 2 . 

In the x « region we seek scaling solutions suggested 
by eq.pl) 



a{x, t) 



®a{x/t a ) 
£l/2-a 



b(x,t) 



<5>b{x/t a ) 

{i/2-a 



(23) 



Several comments are in order to clarify the above scal- 
ing assumptions. First, the scaling of x by the same t a 
in $ a and <!>;, is the assumption that there is only one 
lcngthscale governing the reaction zone. Second, the ex- 
ponent a should be 1/6 or less since the case without 
the membrane gives an upper limit for the spread of the 



reaction zone and there the width is proportional to i 1 / 6 . 
Finally, one should note that the exponent, 1/2 — a, of 
the prefactors of the scaling functions is, in principle, an 
independent exponent. In this case, however, it is fixed 
by the boundary condition (da/dx) x —o ~ and by 

the requirement that the large argument asymptotics of 
b{x, t) should match the solution obtained on the x ~ \ft 
scale. 

Substituting the scaling forms (|2^) into equations 
), one finds that, for large times and for a < 1/2, the 
time derivatives on the left hand sides can be neglected. 
Furthermore, the right-hand sides yield meaningful equa- 
tions only if a is set to a = 1/6. The resulting equations 
then take the form 



(24) 
(25) 





- A$ $ 






dz 2 


n 










dz 2 





where the scaling variable is z = x/t 1 ^ 6 . 

The boundary conditions to the above equations follow 
from a(x — > oo,i) = and db/dx(0,t) = and from 
matching the solutions to the ones found on the diffusive 
scale 



OC =0 , — — 

dz 



2 = 



d$ b 

dz 



z=0 



d$ b 

dz 



b 



(26) 
(27) 



Equations ( |24| . |2 5|) with boundary conditions (|26|,|27|), 

however, pose a difficulty related to the fact that the 
combination v = D a $ a — D^t, satisfies a linear equa- 
tion v" = and the solution, v = Pz + Q, contains a 
integration constant, Q, which is not determined by the 
boundary conditions. Consequently, the scaling functions 
do not appear to be unique. 

This problem of uniqueness can be dealt with by re- 
turning to the diffusive scale, x ~ \/t, and reexamin- 
ing the solutions found there. We shall demonstrate 
the idea on the example of a system where D a = Di, 
(and do = &o since we are at criticality). In this 
case, u = a — b satisfies the diffusion equation for both 
x > and x < 0, the boundary conditions are given by 
it( — oo, t) — —u(oo,t) — ao and d x ii(0~,t) = d x u(0 + ,t) 
and, furthermore, the initial condition [u(x < 0, 0) = ag; 
u(x > 0, 0) = — do] is an odd function of x. It follows then 
that the solution is an odd function, u(x, t) = —u(—x, t). 
Next we note that u = a for x < while u = a — b 
for x > and approaching x = from both sides, the 
oddness of u yields the following relationship 



a{x = Q~,t) = a(x = + , t) - b(x = + , t) 



(28) 



Since there is no accumulation of A particles at x = 0, we 
have da(x = 0~,<) = da(x = + ,i) and, consequently, 
a is continuous function across the membrane, a(x = 
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0~,t) = a{x = 0+,t). Then equation (§|) yields 6(0, t) = 
2a(0, t) which, in turn, provides an additional boundary 
condition for the scaling functions: 

$ 6 (0) = 2$ Q (0) . (29) 

The same extra boundary condition can also be found 
for D a 7^ Db but the argument is rather involved, so 
we shall not reproduce it here. From the perspective of 
critical phenomena, it is quite natural that the scaling 
function does not depend on such details as the diffusion 
coefficients. 

Having the extra boundary condition fl29|), <& Q and $t 
can now be found numerically. Some properties of the 
scaling functions can, however, be seen by just inspect- 
ing the equations. For example, substituting the large 
z asymptotics <&&(z) ~ z into (pi|), one can see that 
<I> Q (z — > oo) is given by the Airy function J2{| . 

On Fig.0, we show that the scaling regime does exist 
and that the numerical results do agree with the solution 
of the full equations (||J|) . 




FIG. 7. Scaling function, <J> a for a (equation (23)). The nu- 
merical solution of the full set of reaction-diffusion equations 
(equations (2) and (3)) is compared with the quasi-stationary 
scaling solution. 

It follows then from (03,E^) that the production rate 
can also be written in a scaling form 

(30) 

and we can observe that the reaction front remains at- 
tached to the wall but it expands with time into the x > 
region. Both the center and the width of the zone diverge 
with time as 

xj ~ w ~ i 1/6 . (31) 

and both exponents are the same in contrast to the de- 
localized phase where Xf ~ t 1 / 2 and w ~ i 1 / 6 . 



The above discussion is based on a mean-field like 
treatment and one can ask what is the role of the fluctu- 
ations which are supposed to be important in low dimen- 
sions. Numerical simulations carried out for an one- and 
two-dimensional systems with semipermeable wall ]30| ] 
indicate that the mean field description discussed above 
remains qualitatively correct. The upper critical dimen- 
sion above which the mean field theory is correct appears 
to be d u — 2. However, in dimension d = 1, the criti- 
cal exponents takes their non mean field values [fl3| . For 
example, the mean position and width exponents at the 
critical point change from the mean field value a = 1/6 
to a = 1/4. 

We can now summarize the properties of the loca- 
lization-delocalization transition discussed above as fol- 
lows. For r < 0, the reaction zone is localized at the 
membrane and the width is determined by the correla- 
tion length, £, describing the penetration of the A par- 
ticles into the constant-density B region. At r = the 
penetration length diverges but there is still a single (di- 
verging with time) length which characterizes the reac- 
tion zone. It should be noted that a diverging diffusion 
length Ijj ~ \/i is always present but it is irrelevant for 
r < 0. For r > 0, however, the diffusion length starts to 
play a role: the reaction zone gets delocalized and two 
distinct lengthscales appear. One of them is the distance 
of the center of the zone from the membrane, xj ~ ^/t, 
which is just the diffusion length while the other is the 
width of the reaction zone, w ~ i 1 / 6 . 

IV. FINAL REMARKS 

The questions of how much C is produced near the 
membrane and whether their density, c, grows to exceed 
some aggregation threshold, cq, may be of importance in 
biological phenomena (e.g. in the building of rather intri- 
cate but regular mineral skeletons of single-cell organisms 
such as radiolaria (27) or diatoms |2q| ). The answers to 
the above questions depend on the localization properties 
of the reaction zone. 

For r < 0, the reaction zone has a finite width and 
thus, provided the C-s do not diffuse away, their density 
will increase with time as c(t) ~ \ft. This result fol- 
lows from the fact that the current, J A (t), of A particles 
towards the reaction zone is proportional to and, 
consequently, the amount of C-s, produced up to time t, 
is given by M c ~ / J A (r)dT ~ y/i. 

A somewhat slower increase of c(t) takes place at r = 0. 
Since the width of the reaction zone diverges asw~ t 1 / 6 , 
one finds c(t) ~ Mc/w ~ i 1 / 3 . We can see that, for both 
r < and r = 0, the density of C-s near the membrane 
exceeds any threshold cq at sufficiently large times. Thus 
supersaturation and, associated with it, the precipitation 
of C may occur in these regimes. 

Finally, for r > 0, the reaction zone leaves the mem- 
brane and only a finite density of C-s left behind. The 
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actual value of this density depends sensitively on the 
initial conditions and one cannot make statements about 
possible precipitation without knowledge of the actual 
parameters. 

The above considerations, of course, do not constitute 
an attempt towards the explanation of a real biologi- 
cal phenomena such as the precipitation of the siliceous 
stuctures of single-cell radiolaria. This is so even if one 
imagines that, at the early stages of the evolution, the 
regular skeletons are either produced as an instability 
in a physico-chemical, reaction-diffusion process or arose 
by surface-tension assisted precipitation where the mem- 
branes are present but play a passive role (their intersec- 
tions defines the precipitation regions) At present 
stage of evolution, the skeletons are covered with a mem- 
branous cytoplasmic sheat which appears to play an im- 
portant role (e.g. transport along the membrane) in the 
skeletal depositions p7j]. Thus any attempt at physico- 
chemical explanation should include the presence of such 
an active membrane near the precipitation zone. 

In this paper, we have derived results for the properties 
of reaction zones near a semipermeable membrane which 
is active only in the sense that it is blocking the trans- 
port of one of the reagents. We hope, however, that our 
results will help discussing more complicated reactions 
near active membranes in the same way as the under- 
standing of the properties of the reaction zone § in the 
A + B — > C reaction helped in elucidating the features 
of the pattern formation in the much more complicated 
Liesegang phenomena 0] . 



ACKNOWLEDGMENTS 

This work has been partially supported by the Swiss 
National Science Foundation in the framework of the Co- 
operation in Science and Research with CEEC/NIS, by 
the Hungarian Academy of Sciences (Grant OTKA T 
019451), and by an EC Network Grant ERB CHRX- 
CT92-0063. Z.R. would like to thank for the hospitality 
of the members of the Theoretical Physics Department 
during his stay at the University of Geneva. 



\-\ 

[8 

[9 
[10 

[11 

[12. 

[13 
[14 

[is; 

[16 

[17; 
[is; 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 
[27 

[28 

[29 

[30 
[31 



H. Taitelbaum, S. Havlin, J. E. Kiefer, B. Trus, and G. 

H. Weiss, J. Stat. Phys. 65, 873 (1991). 

F. Leyvraz and S. Redner, Phys. Rev. Lett. 66, 2168 

(1991); J. Stat. Phys. 65, 1043 (1991). 

E. Ben-Nairn and S. Redner, J. Phys. A25, L575 (1992). 

M. Araujo, S. Havlin. H. Larralde, and H.E. Stanley, 

Phys. Rev. Lett. 68, 1791 (1992). 

H. Larralde, M. Araujo, S. Havlin, and H.E. Stanley. 
Phys. Rev. A46, 855 (1992). 

H. Taitelbaum, Y.E. Lee Koo, S. Havlin, R. Kopelman, 
and GH. Weiss, Phys. Rev. A46, 2151 (1992). 
S. Cornell and M. Droz, Phys. Rev. Lett. 70, 3824 (1993). 
M. Howard and J. Cardy, J. Phys. A28, 3599 (1995). 
GT. Barkema, M.J. Howard and J.L. Cardy, Phys. Rev. 
E53, R2017 (1996). 

S. Cornell, Z. Koza and M. Droz, Phys. Rev. E52, 3500 

(1995) . 

Z. Koza and H. Taitelbaum, Phys.Rev. E54, R1040 

(1996) . 

Z. Koza, J. Stat. Phys. 85, 179 (1996). 

EL. Cabarcos, C-S. Kuo, A. Scala, and R. Bansil, Phys. 

Rev. Lett. 77, 2834 (1996). 

J.S. Langer, Rev. Mod. Phys. 52, 1 (1980). 

GT. Dee, J. Stat. Phys. 39, 705 (1985); Phys.Rev.Lett. 

57, 275 (1986). 

B. Chopard, P. Luthi, and M. Droz, Phys. Rev. Lett. 72, 
1384 (1994); J. Stat. Phys 76, 661 (1994). 
M.C. Cross and PC. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

B. Chopard, M. Droz, and L. Frachebourg, Int. J. Mod. 
Physics, 5, 47, (1994). 

A. Biiki, E. Karpati-Smidroczki, and M. Zrfnyi, J. Chem. 
Phys. 103, 10387 (1995). 

DW. Fawcett, The cell (W.B. Saunders, London, 1981). 
O.R. Anderson, Radiolaria (Springer- Verlag, New York, 
1983). 

Silicon and Siliceous Structures in Biological Systems, 
Eds. T.L. Simpson and B.E. Volcani (Springer- Verlag, 
New York, 1981). 

Handbook of Mathematical Functions, Edited by M. 
Abramowitz and LA. Stegun (Dover, New York, 1965). 
D. Roussin, O. von Suzini, and M. Droz, (unpublished). 
DArcy W. Thompson, On Growth and Form, (Macmil- 
lan, New York, 1942). 



[1] L. Galfi and Z. Racz, Phys. Rev. A38, 3151 (1988). 

[2] Y.E. Lee Koo, L. Li, and R. Kopelman, Mol. Cryst. Liq. 

Cryst. 183, 187 (1990). 
[3] Z. Jiang and C. Ebner, Phys. Rev.A42, 7483 (1990). 
[4] Y-E. Lee Koo and R. Kopelman, J. Stat. Phys. 65, 893 

(1991). 

[5] B. Chopard and M. Droz, Europhys. Lett. 15, 459 (1991). 
[6] S. Cornell, B. Chopard and M. Droz, Phys. Rev. A44 
4826, (1991). 



7 



